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Abstract 



In this paper we study symmetry reductions of a class of nonlinear fourth order partial 
differential equations 

u tt = (ku + iu 2 ) xx + uu xxxx + [iu xxtt + au x u xxx + /3u 2 xx , (1) 

where a, (3, 7, k and [i are arbitrary constants. This equation may be thought of as a 
fourth order analogue of a generalization of the Camassa-Holm equation, about which 
there has been considerable recent interest. Further equation (1) is a "Boussinesq- 
type" equation which arises as a model of vibrations of an anharmonic mass-spring 
chain and admits both "compacton" and conventional solitons. A catalogue of sym- 
metry reductions for equation (1) is obtained using the classical Lie method and the 
nonclassical method due to Bluman and Cole. In particular we obtain several reduc- 
tions using the nonclassical method which are not obtainable through the classical 
method. 



1 Introduction 

In this paper we are concerned with symmetry reductions of the nonlinear fourth order 
partial differential equation given by 

A = u u ~ {ku + iu 2 ) xx - uu xxxx - fJLU xxtt - au x u xxx - (3u 2 xx = 0, (1) 

where a, (5, 7, k and \x are arbitrary constants. This equation may be thought of as an 
alternative to a generalized Camassa-Holm equation (cf. |24| and the references therein) 



u t - eu xxt + 2ku x = uu xxx + auu x + (3u x u xx . (2) 
This is analogous to the Boussinesq equation || [ji]] 

utt = (u xx + \u 2 ) xx (3) 
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which is a soliton equation solvable by inverse scattering jl], 13, 14, 3C, 71], being an 
alternative to the Korteweg-de Vries (KdV) equation 



Ut 



Un 



+ Quu x 



(4) 



another soliton equation, the first to be solved by inverse scattering [39]. 

Two special cases of (|]) have appeared recently in the literature both of which model 
the motion of a dense chain [1621 . The first is obtainable via the transformation 



(u, x, t) I— ► (2ea 3 u + ea 2 , x, t) 
with the appropriate change of parameters, to yield 

utt = {a 2 u + a 3 u 2 ) xx + ea 2 u xxxx + 2ea 3 [uu xxxx + 2u 2 xx + 3u x r 



(5) 



with e > 0. This equation can be thought of as the Boussinesq equation (j||) appended 
with a nonlinear dispersion. It admits both conventional solitons and compact solitons 
often called "compactons" . Compactons are solitary waves with a compact support (cf. 
35|). The compact structures take the form 
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u(x, t) 



3c 2 - 2a 2 



2a 3 



0. 



cos 2 {(12e)- 1 / 2 (x-ci)}, if 

if 



ct\ < 2vr, 
ct\ > 2tt. 



(6) 



or 



u(x, t) 



(7) 



Acos^e)- 1 / 2 x- (|a 3 ) 1/2 i }, if \x - ct\ < 2vr, 
0, if \x - ct\ > 2ir. 

These are "weak" solutions as they do not possess the necessary smoothness at the edges, 



however this would appear not to affect the robustness of a compacton [S2|. Numerical 
experiments seem to show that compactons interact elastically, reemerging with exactly 



the same coherent shape [65]. See [48] for a recent study of non-analytic solutions, in 
particular compacton solutions, of nonlinear wave equations. 

The second equation is obtained from the scaling transformation 

(u,x,t) I— ► (2a 3 u/e, y/e x, t) , 

again with appropriate parameterisation, 

u t t = {a 2 u + ot3U 2 ) xx + eu xxtt + 2ea 3 [ 



+ 2u xx + 3u x u x 



(8) 



with e > 0. This equation, unlike (|5|) is well posed. It also admits conventional solitons 
and allows compactons like 



u(x, t) 



Ac 2 - 3a 2 



2a 3 



cos 



0. 



|(12e)" 1/2 (x -ct)\ , if \x-ct\<2TT, 

if \x — ct\ > 2"7T, 



(9) 



or 



u(x, t) 




\l/2 . 
a 2 ) t 



| , if \x — ct\ < 2tt, 
if \x — ct\ > 2ir. 



(10) 
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These again are weak solutions, and are very similar to the previous solutions: both (Q) 
and (|l^) are solutions with a variable speed linked to the amplitude of the wave, whereas 
both (|6|) and (||) are solutions with arbitrary amplitudes, whilst the wave speed is fixed 
by the parameters of the equation. 

The Fuchssteiner-Fokas-Camassa-Holm (FFCH) equation 

ut — u xx t + 2ku x = uu xxx — 3uu x + 2u x u xx , (11) 

first arose in the work of Fuchssteiner and Fokas [p4| , |36| ] using a bi-Hamiltonian approach; 
we remark that it is only implicitly written in - see equations (26e) and (30) in this 



paper - though is explicitly written down in [34]. It has recently been rederived by Camassa 



and Holm [11] from physical considerations as a model for dispersive shallow water waves. 



In the case k = 0, it admits an unusual solitary wave solution 
u(x, t) = Aexp (—\x — ct\) , 



where A and c are arbitrary constants, which is called a "peakon". A Lax-pair [11] and bi- 



Hamiltonian structure [36 have been found for the FFCH equation (|Tl]) and so it appears 
to be completely integrable. Recently the FFCH equation (11) has attracted considerable 



attention. In addition to the aforementioned, other studies include [12, 25, 26, 27, 29, 32 



ii m m mm- 



Symmetry reductions and exact solutions have several different important applications 
in the context of differential equations. Since solutions of partial differential equations 
asymptotically tend to solutions of lower-dimensional equations obtained by symmetry 
reduction, some of these special solutions will illustrate important physical phenomena. 
In particular, exact solutions arising from symmetry methods can often be used effectively 
to study properties such as asymptotics and "blow-up" (cf. |3^, [38|). Furthermore, explicit 
solutions (such as those found by symmetry methods) can play an important role in the 
design and testing of numerical integrators; these solutions provide an important practical 
check on the accuracy and reliability of such integrators (cf. [||, |67|). 

The classical method for finding symmetry reductions of partial differential equations 
is the Lie group method of infinitesimal transformations, which in practice is a two-step 
procedure (see § 2 for more details). The first step is entirely algorithmic, though often 
both tedious and virtually unmanageable manually. As a result, symbolic manipulation 
(SM) programs have been developed to aid the calculations; an excellent survey of the 
different packages available and a description of their strengths and applications is given 



by Hereman [41]. In this paper we use the MACSYMA package symmgrp.max [15] to calcu- 
late the determining equations. The second step involves heuristic integration procedures 
which have been implemented in some SM programs and are largely successful, though not 
infallible. Commonly, the overdetermined systems to be solved are simple, and heuristic 
integration is both fast and effective. However, there are occasions where heuristics can 
break down (cf. |5(J for further details and examples). Of particular importance to this 
study, is if the classical method is applied to a partial differential equation which con- 
tains arbitrary parameters, such as (||) or more generally, arbitrary functions. Heuristics 
usually yield the general solution yet miss those special cases of the parameters and ar- 
bitrary functions where additional symmetries lie. In contrast the method of differential 
Grobner bases (DGBs), which we describe below, has proved effective in coping with such 
difficulties (cf. @, §(| ||). 
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In recent years the nonclassical method due to Bluman and Cole [0] (in the sequel 
referred to as the "nonclassical method" ) , sometimes referred to as the "method of partial 
symmetries of the first type" [38|, or the "method of conditional symmetries" |45], and the 
direct method due to Clarkson and Kruskal fHj^ have been used, with much success, to gen- 
erate many new symmetry reductions and exact solutions for several physically significant 
partial differential equations that are not obtainable using the classical Lie method (cf. 
fli~6| , |i~9f | and the references therein). The nonclassical method is a generalization of the clas- 
sical Lie method, whereas the direct method is an ansatz-based approach which involves 
no group theoretic techniques. Nucci and Clarkson [53] showed that for the Fitzhugh- 
Nagumo equation the nonclassical method is more general than the direct method, since 
they demonstrated the existence of a solution of the Fitzhugh-Nagumo equation, obtain- 
able using the nonclassical method but not using the direct method. Subsequently Olver 



5[ (see also || 57]) has proved the general result that for a scalar equation, every reduc- 
tion obtainable using the direct method is also obtainable using the nonclassical method. 
Consequently we use the nonclassical method in this paper rather than the direct method. 

The method used to find solutions of the determining equations in both the classical 
and nonclassical method is that of DGBs, defined to be a basis 13 of the differential ideal 
generated by the system such that every member of the ideal pseudo-reduces to zero with 
respect to the basis 6. This method provides a systematic framework for finding inte- 
grability and compatibility conditions of an overdetermined system of partial differential 
equations. It avoids the problems of infinite loops in reduction processes and yields, as far 
as is currently possible, a "triangulation" of the system from which the solution set can be 



derived more easily 52, 61]. In a sense, a DGB provides the maximum amount of 



information possible using elementary differential and algebraic processes in finite time. 

In pseudo-reduction, one must, if necessary, multiply the expression being reduced 
by differential (non-constant) coefficients of the highest derivative terms of the reducing 
equation, so that the algorithms used will terminate |52|]. In practice, such coefficients 
are assumed to be non-zero, and one needs to deal with the possibility of them being zero 
separately. These are called singular cases. 

The triangulations of the systems of determining equations for infinitesimals arising 
in the nonclassical method in this paper were all performed using the MAPLE package 
diffgrob2 pi. This package was written specifically to handle nonlinear equations of 
polynomial type. All calculations are strictly 'polynomial', that is, there is no division. 
Implemented there are the Kolchin-Ritt algorithm using pseudo-reduction instead of re- 
duction, and extra algorithms needed to calculate a DGB (as far as possible using the 
current theory), for those cases where the Kolchin-Ritt algorithm is not sufficient [52]. 
The package was designed to be used interactively as well as algorithmically, and much 
use is made of this fact here. It has proved useful for solving many fully nonlinear systems 



(cf. 0,0|g,p|2|). 

In the following sections we shall consider the cases /i = and /i ^ 0, when we set fi = 1 
without loss of generality, separately because the presence or lack of the corresponding 
fourth order term is significant. In § 2 we find the classical Lie group of symmetries and 
associated reductions of ([]]). In § 3 we discuss the nonclassical symmetries and reductions 
of (I) in the generic case. In § 4 we consider special cases of the the nonclassical method 
in the so-called r = case; in full generality this case is somewhat intractable. In § 5 we 
discuss our results. 
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2 Classical Symmetries 

To apply the classical method we consider the one-parameter Lie group of infinitesimal 
transformations in (x, t, u) given by 

x * =x + e£ (x, t,u) + (e 2 ) , 

t* = t + er(x,t,u) +0 (e 2 ) , (12) 

u* = u + e(j){x, t,u)+0 (e 2 ) , 

where e is the group parameter. Then one requires that this transformation leaves invariant 
the set 

S A = {u(x,t) : A = 0} (13) 

of solutions of (|l|). This yields an overdetermined, linear system of equations for the 
infinitesimals £(x,t, u), r(x,t,u) and 4>(x,t,u). The associated Lie algebra is realised by 
vector fields of the form 

3 d 

Having determined the infinitesimals, the symmetry variables are found by solving the 
characteristic equation 

dx _ dt _ du 
£(x,t,u) r(x,t,u) (f>(x, t, it)' 

which is equivalent to solving the invariant surface condition 

ip = t, u)u x + r(x, t, u)u t — (p(x, t, u) = 0. (16) 
The set Sa is invariant under the transformation ((TJ) provided that 

pr( 4 )v(A)| As o = 0, 

where pr^ 4 W is the fourth prolongation of the vector field (|l4|), which is given explicitly 
in terms of £, r and <p (cf- [54]). This procedure yields the determining equations. There 
are two cases to consider, (i) \i = and (ii) /j / 0. 

2.1 Case (i) fi = 0. 

In this case we generate 15 determining equations, using the MACSYMA package 
symmgrp .max. 

t u = 0, t x = 0, £ u = 0, 4> uu = 0, & = 0, 
a(4> u u - (j>) = 0, fi{(f> u u - 4>) = 0, 2<j) tu - T tt = 0, 
tyxuU - 6£ xx u + a</> x , 2r t u - 4£ x u + = 0, 
4P<frxu + 3a4> xu ~ 2fi£,xx ~ 3a£ xx = 0, 

4>tt - <t>xxxxU ~ 2j4>xxU ~ K<p xx = 0, 

3a(f) xxu u + 2^(f) u u + 4£c7u - a^ xxx u - 2-ycj) = 0, 

6<f>xxuU 2 + 4^7« 2 - 4£ xxx u 2 + 2(3(f) xx u + 2£ x ku - K<f> = 0, 

HxxxuU + 4j<p xu u - 2£ xx ju - i xxxx u + a(f) xxx + Aj(t) x + 2k(\> xu - ^ xx k = 0, 
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and then use reduceall in dif f grob2 to simplify them to the following system 

Zu = 0, 6 = 0, 7(14/9 + 90)^ = 0, r u = 0, 
7/<14/3 + 9a)r t = 0, = 0, 7re(14/3 + 9a)0 = 0. 

Thus we have special cases when 7 = 0, re = and/or 14/3 + 9a = 0. The latter condition 
provides nothing different unless we specialize further and consider the special case when 
a = — I and (3 = ||. We continue to use reduceall in diffgrob2 for the various 
combinations and it transpires that there are only four combinations which yield different 
infinitesimals. Where a parameter is not included it is presumed to be arbitrary 

(a) re = 0, 

(b) 7 = 0, 

( C ) 7 = K = 0, 



(d) a 



5 
2' 



/3=§, 7 = re = 0, 





= 0, 


6 


= 0, 


tx — 0, T u — 


0, 






= 0, 




= 0, 


2r t u + <p = 0. 






£« 


= 0, 




= 0, 


£,xx = 0, T u = 


= 0, 




6, 


-T t = 


0, 




= 0, 2&.u + 


= 0. 




& 


= 0, 




= 0, 


Cxx = 0, T u = 


= 0, 




Ttt 


= 0, 




= 0, 


2r t ii - 4£ x u - 




0. 


tu 


= 0, 




= 0, 


€xxx = 0, T u 


= 0, 




Ttt 


= 0, 




= 0, 


2r t u - 4£ x u - 


= 


0. 



Hence we obtain the following infinitesimals. 

Table 2.1 



Parameters 


e 


T 


4> 






Cl 


C2 





(2-li) 


K = 


Cl 


C3t + C2 


— 2C3« 


(2.1ii) 


7 = 


C3X + Cl 


c%t + C 2 


2C3U 


(2.1iii) 


7 = re = 


C4X + Cl 


+ C 2 


(4c 4 - 2c 3 )u 


(2.1iv) 


Ui 2, l~> 28 

7 = re = 


C5X 2 + C4X + Cl 


+ C 2 


[4(2c 5 x + c 4 ) - 2c 3 ]u 


(2.1v) 



where c\, C2, . . . ,05 are arbitrary constants. 

Solving the invariant surface condition yields the following seven different canonical 
reductions: 

Reduction 2.1 a, 8, 7 and re arbitrary. If in (2. li-2.lv) £3 = 04 = 05= 0, then we may 
set C2 = 1 without loss of generality. Thus we obtain the reduction 



u(x, t) = w(z), 
where w(z) satisfies 



(re - of) 



d 2 w 
dl2 



+ 2 7 



Z = X — C\t 



d 2 w ( dw \ 2 

, ^z 2 + \~dz~) 



d 4 w dw d 3 w ( d 2 w 
+ w—r-r + a— ;— 3- + 8 



dz 4 



dz dz 3 



dz 2 



= 0. 
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Reduction 2.2 a, (3 and 7 arbitrary, k = 0. If in (2.1ii), (2.1iv) and (2.1v) 04 = 05 = 0, 
C3 / 0, then we may set C2 = and C3 = 1 without loss of generality. Thus we obtain the 
reduction 



u(x, t) = t 2 w(z), 
where w(z) satisfies 



z = x-ci log(i), 



d 4 w dw d 3 



w 



+ a 



w 



dz 4 
+2 7 



dz dz 3 
d 2 w 



+ P 



d 2 w 
~dz* 



w 



dz 2 



+ 



dw 
dl 



d 2 w dw 

TT ~ 5c i j Gw 

dz z dz 



0. 



Reduction 2.3 a, [5 and k arbitrary, 7 = 0. If in (2.1iii) C3 / 0, then we may set 
ci = C2 = and C3 = 1 without loss of generality. Thus we obtain the reduction 



u(x, t) = t 2 w(z), 
where w{z) satisfies 



x/t, 



d 4 w dw d 3 w ( d 2 w 
w—r-r + ct— 5—5- + p 



dz 4 



dz dz 3 



dz 2 



d 2 



w 



,d 2 



w 



dw 



dz 2 



. 9 + 2z— -2w = 0. 

dz z dz 



Reduction 2.4 a and /3 arbitrary, k = 7 = 0. If in (2.1iv) and (2.1v) C3 = C5 = and 
C4 / 0, then we may set c\ = and C2 = 1 without loss of generality. Thus we obtain the 
reduction 



u(x,t) = w(z) exp(4c4^), 
where w(z) satisfies 



z = x exp(— C4t), 



d 4 w dw d 3 w 



w 



dz 4 



+ a 



dz dz 3 



+ P 



f d 2 w 
[jz 2 



2 od 2 W 2 dw 2 

CaZ -—r + icaZ- 16c 4 w = 0. 

dz z dz 



Reduction 2.5 a and f3 arbitrary, k = 7 = 0. If in (2.1iv) and (2.1v) C5 = and C3C4 / 0, 
then we may set c\ = C2 = and C3 = 1 without loss of generality. Thus we obtain the 
reduction 



u(x,t) =w(z)t 4c4 ~ 2 , 
where w{z) satisfies 



z = xt~ c \ 



dw d 3 w 



d 4 „ 
w—r-r + a 



dz 4 dz dz 3 
,d 2 u> 



w 



dz 2 



clz 2 %^r + (7c 2 . - 5c 4 ) z^- - (16c 2 . - 20c 4 + 6) w = 0. 

dz z ' dz ' 
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Reduction 2.6 a = — |, (3 = ||, 7 = k = 0. If in (2.1v) C3 = and C5 / 0, then we may 
set ci = — mc, C2 = 1, C4 = and C5 = m/c, without loss of generality. Thus we obtain 
the reduction 

w(z) exp(— 8mt) fx — c\ 
u(x,t) = — — — — ^, z = exp(-2mt), 



[z - exp(-2mt)] 8 \x + c 

where w(z) satisfies 

d 4 w dw d 3 w f d 2 u> x 2 
28w—r - 70- r^r + 45 



dz 4 dz dz 3 \ dz 

-c 4 m 2 ( 1792z 2 ^ - 12544z^ + 28672uA = 0. 
\ dz z dz J 

Reduction 2.7 a = — §, j3 = ||, 7 = k = 0. If in (2.1v) C3C5 / 0, then we may set 
ci = — mc, C2 = C4 = 0, C3 = 1 and C5 = m/c, without loss of generality. Thus we obtain 
the reduction 

. . w(z)t- 2 ^ 1+Am ) (x-c\, m 
u(x,t) - - - 1 ' • -'• 



(z-t~ 2m f ' \x + c 

where w(z) satisfies 

„„ d 4 w „ A dw d 3 w A „ f d 2 w\ 2 4 2 2<i 2 w 

28w— r - 70- + 45 — r - 1792c 4 m V 



dz 4 dz dz 3 \ dz 2 J dz 2 

dw 
~dz~ 



+ (12544m 2 - 4480m) c 4 z^- - (28672m 2 - 17920m + 2688) c 4 u> = 0. 



2.2 Case (ii) fj, ^ 0. 

In this case we generate 18 determining equations, 

T„ = 0, T x = 0, £ M = 0, & = 0, (p uu = 0, 

xtu = 0, a{4> u u - <f>) = 0, 20 tu - r tt = 0, /3(0„u - 0) = 0, 

20a;« - ixx = 0, 40 xu u - 6£ xx u + a(f> x = 0, 

2r t u - 2£ x u + = 0, 

4/30a;« + 3a0 xtt - 2/3£ xx - 3a£ xx = 0, 

<^xxtt« + 2r t u - 4£ x u + = 0, 

3a4> xxu u + 270 u u + 4^7« - a£ xxx u - 2<ycf) = 0, 

<Ptt ~ 4>xxxxU - 2j4> xx u - K<p xx - 4> xxtt = 0, 

64>xxuU 2 + 4^7u 2 - A£ xxx u 2 + 2(3<t> xx u + ttu u + 2£ x fiu - K = 0, 

40 :ra: r U U + ^4>xuU - 2£ xx ju - £ XXXX U + CX(f> xxx + 470^ + 2n<j) xu - £ XX K = 0, 
and then use reduceall in dif f grob2 to simplify them to the following system, 

Cu = 0, & = 0, £ x = 0, T M = 0, KT t = 0, T x = 0, K0 = 0. 
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Here k = is the only special case, yielding the slightly different system 

£« = 0, 6 = 0, £s = 0, Tu = 0, T tt = 0, t x = 0, + 2r t u = O. 

Thus we have two different sets of infinitesimals, and in both cases a, (3 and 7 remain 
arbitrary. 



Table 2.2 



Parameters £ 


r 







Cl 







(2.2i) 


K = Cl 


C3t + C 2 


-2c 3 u 


(2.2ii) 



From these we have the following two canonical reductions: 

Reduction 2.8 a, 0, 7 and k arbitrary. If in (2.2i) and (2.2ii) C3 = 0, then we may set 
C2 = 1 without loss of generality. Thus we obtain the following reduction 



u(x, t) = w(z), 
where w(z) satisfies 



(« - 4) 



d 2 w 

dl2 



+ 27 



X — C\t, 



d 2 w ^ ( dw x 
dz 2 \ dz 



d A w dw d 3 w ( d 2 w 
+w- rir + a- — — + p 



dz 4 



dz dz 3 



dz 2 



+ 4 



d A w 



0. 



Reduction 2.9 a, f3, 7 arbitrary, k = 0. If in (2.2ii) C3 7^ 0, then we may set C2 = 0, 
C3 = 1 without loss of generality. Thus we obtain the following reduction 



u(x, t) = t 2 w(z), 
where w(z) satisfies 



x - ci log(i), 



d 4 w dw d 3 w ( d 2 w 
w—r^r + a— ;— 5- + p 



dz 4 



dz dz 3 



dz 2 



+ 4 



d 4 



dz 4 



d 3 w 
+5ci^3+2 7 



dw 



d 2 w 

dz 2 V dz 



+ (6 - c 2 ] 



d 2 



w 



dz 2 



dw 

5ci— 6w = 0. 

dz 



2.3 Travelling wave reductions 

As was seen in § 1, special cases of (|]) admit interesting travelling wave solutions, namely 
compactons. In this subsection we look for such solitary waves and others, in the frame- 
work of (|l|). Starting with compacton-type solutions, we seek solutions of the form 



u(x, t) = a 2 cos n {a3(x — ait)} + 04, 



(17) 
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where 01,02,03,04 are constants to be determined. We include the (possibly non-zero) 
constant 04 since u is open to translation. The specific form of the translation will put 
conditions on 04, which may or may not put further conditions on the other parameters 
111 (!U) and those in (@) - see below. If n = 1 we have the solutions, where the absence of 
a parameter implies it is arbitrary, 



(i) a 



111 



a 



0, 
1, 

a 



-1, 7 = 0, a 4 



K — OT 



2 2 



f3 = 0, 7 >0, a 2 (l + 2 7 //) 
1/0 



- > 0, al 

a 6 



27 

■ 

a 



04 



K = 0, a 3 = 27. 
o(k — a 2 ) — 2a^7/i 



2 7 (1 - a) 

These become n = 2 solutions via the trigonometric identity cos 26 = 2 cos 2 6 — 1. By 
earlier reasoning the associated compactons are weak solutions of (||). When considering 
more general n we restrict n to be either 3 or > 4 else the fourth derivatives of u(x, t) that 
we require in (Q) would have singularities at the edges of the humps; we find 

2^2 — n n, . 97 

-, 0= , 7 >0, a?(l + 2 7 /x)-K = 0, al = ±, 

n 



a 



04 



-/UOi- 



n n 
When n = 3 or n = 4 our compacton would be a weak solution since not all the derivatives 
of u(x, t) in (H) in these instances are continuous at the edges. For n > 4 the solutions are 
strong. 

For more usual solitary waves we seek solutions of the form 



u(x,t) = a2sech n {a3(x — ait)} + 04, 

where a±, 02, 03, 04 are constants to be determined. If n = 2 then a = —1, (3 
have solutions 

(i) 7 <0, o?(l + 2 7M )- 



-2 and we 



0, 



57, a 4 = -jr(2a 2 + 3ai//), 



(ii) 7 2 /4a|, a 2 



3a|(« 



2of7/z) 



(2a|- 7 )(2a| + 7) 
and for general n, including n = 2 (7 > 0) 
2 „ n+2 2 



04 



2(2a|+ 7 ) 



a 







ai(l + 27M) 



0, 



«3 



7 



04 



n n 

Now consider the general travelling wave reduction, u(x, t) = w(z), z = x — ct, where w(z) 
satisfies 



d 2 w 
dl 2 



+ 2 7 



d 2 w 
'dz 2 ' 



( dw x 
V d7 



, d 4 u> 



d 4 w; dw d 3 u; 
+ ^-r-r + Q!— r^r + P 



d 2 wY 

dz 2 J 



0. 



dz 4 dz 4 dz dz 3 
In the special case (3 = a — 1, we can integrate this twice with respect to 2 to give 



/ 2\ 2 2 ^ ^ ^ / 

[K - c j w + + fic -—5- + -(a 

dz z 2 



— + Az + B 
dz J 







(18) 
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with A and B the constants of integration. If we assume A = 0, then we make the 



transformation W(z) = w(z) + fic 2 , multiply (18) by W 
to z to yield 



a-3 



dW 



and integrate with respect 



l W a + -^-w a - 1 + 



A 



a 



a — 1 



2 w a ~ 2 + W a ~ 2 



a 



dW 
~dz~ 



+ C = 



for q / 0,1, 2, where Ai = k — c 2 — 27/xc 2 and A2 = 1? — ^c 2 (k — c 2 — 7/ic 2 ). In the 
special cases a = 0, 1, 2 we obtain respectively 



. __. Ax 2A 2 1 /diy 
7 lo g ^- — -— T + — 5 



W VF 2 W 2 \ dz 
+ A x log 



A 2 1 fdW 



W + W V dz 



1 n /dW 
-7^ 2 + AxW + A 2 logW+ — 

2 V dz 



+ C = 0, 



+ C = 0, 



+ C = 0, 



(19) 
(20) 
(21) 



where C is a constant of integration. For a £ Z, an integer, with a > 3, (|2.3[) may be 
written as 



(w + 



2 nQ-2 / duA 



dz J 



+ 



+- 



jw 
a 

a — 



2\ a - 



(w + fjLC ) 



2 a (k — c 2 ) — 2/L17C 2 



.2 a ~ 3 f,.JZ\a-3-n 



+ 



£ 

n=0 



(n + 2) 



Q: 



(22) 



2\« — 2 „ 

+ /ic 2 ) + C 



(/ic 2 ) " 1 [a (as - c 2 ) - 2^ 7 c 2 
a(a — l)(a — 2) 



0. 



where C is a constant of integration. If we require that it; and its derivatives tend to zero 
as z — > ±00, then B = D = 0. Ifa = 3 this equation induces so-called peakons (cf. p| ) 
as a (/t - c 2 ) — 2/i7c 2 — ► (see [12, |(], 44, |6^|). Similarly if a = 4 this equation is of 
the form found in [40| which induces the 'wave of greatest height' found in [31|. Both 
solutions, in their limit, have a discontinuity in their first derivative at its peak. Note that 
if q (k — c 2 ) — 2^7c 2 = 0, equation ( [22] ) becomes 



(w + \xc 



,a-2 



— + ^ 2 
dz J a 



0. 



(23) 



Since a > then we require 7 < to give a peakon of the form 



u(x, t) 



a c 



k) 



27 



exp 



a 



1/2 



\x — ct\ 



(24) 



The height of the wave, because of the form of (|l]), is dependent upon the square of the 
speed, whereas the peakons in [11] and [31] are proportional to the wave speed. 
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3 Nonclassical symmetries (r ^ 0) 

In the nonclassical method one requires only the subset of Sa given by 

Sa,^ = {u{x, t) : A(u) = 0, i/j(u) = 0}, (25) 
where Sa is defined in ( |l3| ) and ifi = is the invariant surface condition (16), to be 



invariant under the transformation ([12]). The usual method of applying the nonclassical 
method (e.g. as described in f45|), involves applying the prolongation pi^v to the system 
composed of (|l|) and the invariant surface condition ( |l6D and requiring that the resulting 
expressions vanish for u 6 Sa,-^, i.e. 

pr( 4 )v(A)| {A=0il/;=0} = 0, pr (1) v(V)| { A=o^=o} = 0. (26) 

It is well known that the latter vanishes identically when ip = without imposing any 
conditions upon £, r and (p. To apply the method in practice we advocate the algorithm 
described in p2| for calculating the determining equations, which avoids difficulties arising 



from using differential consequences of the invariant surface condition (16). 



In the canonical case when r^Owe set r = 1 without loss of generality. We proceed 
by eliminating utt and u xx tt in (|l]) using the invariant surface condition (16) which yields 

££xU x + 2u 2 ^ u - 2(p u £u x + £ 2 u xx - (f> x £ - i t u x + (fxpu ~ 4>£,uU x + <pt~ ku xx 



2(p xx Cx - 2(p xu (j) x - 4££v,. 



X^XX 



-2j(uu xx + u 2 x ) - uu xxxx - au x u xxx - (3u 2 xx + fi 

^Ptu^xx tpu^Pxx <Pxxt £,xxtUx 4 > tuuU x ~\~ 4>x£,xx ~\~ CtuuU x ~\~ ^t^xxx £ "xxxx 
~\~<Pxxx£, ^Pu^xx (pfixxu CCxxx^x 'PCuuu'U'x 'PCxxuUx ~\~ (p^u^xxx 4 ) 4 ) uu'Uxx 

4 > 4 > uuuU x ~\~ 2£ >x tU xx + 2t^ x t u U x 2(p x t u ll x Q£,x£,xx'U'X ^Cu£,xxU x ^Cxx^xx 
~^~2(p u ^ xx ll x §£,uu£,x'U'x &£,xuCx'U'x ^^>CuCxUxU xx <l£,CxUxxx ~i~ ^^Pu^x^xx (27) 
~\~^4 > uu£,x'U' X ~\~ 6(pxu£,xUx 2^ UUU U X 5£,£,xuuU x ~\~ 2(p^ xuu U x ^Cu^,uuU x 

12£ I !^ UU U X U XX + ^(p^uuU x U xx + 4<^> n £ mi - U 2 . -|- ^(px^uu^x ^^Cxxu^x ^^Cu^xu^x 

^^CCxu^x^xx 2(p£ >xu 11 xx + 6 < Pu£,xu'U'x ^(frxCxullx 12£ iU U x U xx 8^ U U X U XXX 

-6^ u ul x + 9<p u £ u u x u xx + 3<p x ^ u u xx + 5(p uu £uul + 8<Pxu£uul + 3<p xx £ u u x 

~ > <~'^^tuU x U xx 

XXX x^xx + 5><Pxuiu xx + 2<p 

ami 

£u x + 5(p X uu£,U 2 



X 



'xxuC^x 3(p u (p uu U x 2(p uu <p x U x 2<p(p xuu ll x 4:(pu4 > xu'U'X 



0. 



We note that this equation now involves the infinitesimals ^ and <p> that are to be deter- 
mined. Then we apply the classical Lie algorithm to ( p7| ) using the fourth prolongation 
p r ( 4 ) v and eliminating u xxxx using (p7|). It should be noted that the coefficient of u xxxx 
is (£ 2 + ixu). Therefore, if this is zero the removal of u xxxx using ( |27| ) is invalid and so 
the next highest derivative term, u xxx , should be used instead. We note again that this 
has a coefficient that may be zero so that in the case \i / and i 2 + fiu = one needs 
to calculate the determining equations for the subcases non-zero separately. Continuing 
in this fashion, there is a cascade of subcases to be considered. In the remainder of this 
section, we consider these subcases in turn. First, however, we discuss the case given 
by m = 0. 
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3.1 Case (i) fi = 0. 

In this case we generate the following 12 determining equations. 

e« = o, 

^xuuuU + 3ac/) xuu = 0, 
6<t>uuuU + 2f3(f) uu + 3a(p uu = 0, 

^4>uuU 2 + Ot(j) u U — Ct(j) = 0, 

4(f> xu u - 6£ xx u + a<f> x = 0, 
3<f>uuU 2 + (3(f) u u - /3(f) = 0, 

120 xmt it + 4/?0 xu + Sacf) xu - 2(3£ xx - 3a£ xx = 0, 

uuU 

+3acf) xxu u + 2-y(f) u u + 4^7u - a£ xxx u - 2-/(f> = 0, 
6<Pxxuii 2 + 4^ 7 u 2 - 4^ x:r:r u 2 + 2/3(f) xx u + 2&.KU - 4£ 2 £ r « - 2& t u - K<f> + £ 2 ^ = 0, 

- 4>xxxxU 2 - 2j(p xx u 2 - K(f) xx U - 4££ X (/) X U - 2£ t <fixU 

+4> 2 <PuuU + 4£ E 00 u u + 2(p(p tu u + 4£ x </> t u + £# x - 4> 2 <Pu ~ <p(j>t = 0, 

^4>xxxuU 2 + 4j(/) xu U 2 - 2£ xx ~fU 2 - £ XXXX U 2 

+a(f) xxx u + 4j(/> x u + 2£(f><f> uu u + 2n<j) xu u + 8^ x <p u u + 2£ t ^ M -u 
+2£<f) tu u - £ xx ku - 4££ 2 u + 2^ t ^n + - 2^^> u + ££ x (f> - £ t (f) = 0. 

As guaranteed by the nonclassical method, we get all the classical reductions, but we also 
have some infinitesimals that lead to nonclassical reductions, namely 



Table 3.1 


Parameters 








K = 

a = (3 = 7 = 





d 2 g d<7 o 

q(t)u where —— r + o- q =0 

yv ; dt 2 y di 

C32/ 3 + c 2 y 2 + ciy + c (y = x±- v /Kt) 

7-T-— + (C4X +C32 +C2X +C1X + C0J 

flf(t) dt 


(3-li) 
(3-lii) 


a = /3 = 7 = K = 





where 

2 d 3 g d 5 d 2 a /d 5 \ 3 4 

5 d^- 4ff d^ +2 U] +24C49 =0 


(3.1iii) 



From these we obtain three canonical reductions. 

Reduction 3.1 a, (3, 7 arbitrary, k = 0. In (3. li) we solve the equation for g(t) by 
writing g(t) = [ln(tp(t))]t then ip(t) satisfies 

(^) 2 = 4c^ 3 + c 2 , (28) 
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where c\ and C2 are arbitrary constants; c\ = C2 = is not allowed since g(t) ^ 0. Hence 
we obtain the following reduction 

u(x, t) = w(x)ip(t), 

where w(x) satisfies 

^4„„ J„„ j3 



d A w ^ dw d^w ^ ( d 2 w \ 2 ^ ^ 
dx 4 dx dx 3 \ dx 2 / 



w 



d 2 w ( dw \ 2 
+ — 



6c\w = 0. 



dx 2 \ dx / 

There are three cases to consider in the solution of (|2q) . 

(i) If ci = 0, we may assume that ip(t) = t without loss of generality. 

(ii) If C2 = 0, then ip = [c2(t + C3) 2 ] 1 and we may set C2 = 1, C3 = without loss of 
generality. 

(iii) If c\C2 7^ we may set c% = 1, C2 = —53 without loss of generality so that if)(t) is 
any solution of the Weierstrass elliptic function equation 



dp 
dt 



(t;0, 93 ) 



4p 3 (t;0, 53 ) -53- 



(29) 



Reduction 3.2 k arbitrary, a = (3 = 7 = 0. From (3.1ii) we get the following reduction 



u(x,t) = w(z) ± — -^y 4 ± ^=y 3 ± + 0)*, y = x±V^t, z = x=F^t, 

8a/ k 6a/k 4a/k 

where w{z) satisfies 

,— d 4 w 
^dx^ ±3C3 = °- 
This gives us the exact solution 

u(x, t) = =f-^=z 4 + c A z 3 + c 5 z 2 + cqz + c 7 ± ^=y 4 ± 7-7=y 3 ± T-7^y 2 + c o*- 

Oa/k oa/k Da/k 4a/k 

Reduction 3.3 a = /? = 7 = K = 0. In (3.1iii) we integrate our equation for g{t) up to 
an expression with quadratures 

2 



C2 „ .3 C l „,2 



d 2 / d \ 
5 d^" 2 (df) + 24c W 5 2 (s)ds + 24c 55 = 0. 



(30) 



We get the following reduction 



u(x, t) 



1 



9(t) 

where w(x) satisfies 
d A w 



w(x) + (C4X 4 + C3X 3 + C2X 2 + C\X + Cq 



) f9\s) 



ds 



dx 4 



24cs = 0. 
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This is easily solved to give the solution 
1 



u(x, t) 



C 5 X + CqX 3 + c 7 x 2 + CgX + c 9 



9(t) 

+ (c 4 x 4 + c 3 x 3 + c 2 x 2 + c\x + c ) J g 2 (s) ds 
where g(t) satisfies (|30"1). 



3.2 Case (ii) fj, ^ 0. 

As discussed earlier in this section, we must consider, in addition to the general case of 
the determining equations, each of the singular cases of the determining equations. 

3.2.1 £ 2 + u ^0. 

In this the generic case we generate 12 determining equations - see appendix A for de- 
tails. As expected we have all the classical reductions, however we also have the following 
infinitesimals that lead to genuine nonclassical reductions (i.e. not a classical reduction). 

Table 3.2 



Parameters 









K = 





d 2 g dg 
g[t)u where + g— -g = 


(3.2i) 


1 + 2 7 = 


C\t + C2 


-2c\(c\t + C2) 


(3.2ii) 


K = 1 + 27 = 


c 2 (t + ci) 2 


u(t + C1)- 1 -3c|(t + ci) 3 


(3.2iii) 


a = (3 = 7 = 




C32/ 3 + c 2 y 2 + ciy + c (y = x ± y/nt) 


(3.2iv) 


a = -§,/3 = 2,7 = 


±\y/K(x + Cl) 


±2^u± i/€ 3 / 2 (x + ci) 2 


(3.2v) 


a = (3 = j = K = 





c 3 x 3 + c 2 x 2 + ClX + c 


(3.2vi) 


a = (3 = j = K = 





(u + C3X 3 + C2X 2 + C\X + Co)(t + C4) _1 


(3.2vii) 



From these infinitesimals we obtain six reductions. 

Reduction 3.4 a, (3, 7 arbitrary, k = 0. In (3.2i) we solve the equation for g(t) by 
writing g(t) = [ln(tp(t))]t then tp(t) satisfies 



(^) 2 = 4c^ 3 + c 2 (31) 

though ci = C2 = is not allowed to preserve the fact that g(t) ^ 0. We obtain the 
following reduction 

u(x, t) = w(x)ip(t), 
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where w(x) satisfies 



2„„\ 2 



d 4 w dw d 3 w „ ( d 2 w . 
ax dx 6 \ dx z 



dx 4 



d 2 w f dw\^ 
dx 2 V dx / 



, d 2 w 

!6c l( ^ -.]-.). 



There are three cases to consider in the solution of (|3] 

(i) If ci = 0, we may assume that ip(t) = t without loss of generality. 

(ii) If C2 = 0, then -0 = \p%(p + C3) 2 ] 1 and we may set C2 = 1 and C3 = without loss 
of generality. 

(iii) If c\C2 7^ we may set ci = 1 and C2 = —53 without loss of generality so that ip(t) 
is any solution of the Weierstrass elliptic function equation p3). 

Note that in the special case 
d 2 w 

w = 0, 



dz 2 



we are able to lift the restrictions on ip(t) so that it is arbitrary, if (3 + 1 + 27 = a + 27 = 0. 
This yields the exact solution 

u(x,t) = ip{t) (c 2 e x + c 3 e~ x ) , 

where ip(t) is arbitrary, k = 0, a = —27 and (3 = —1 — 27. 

Reduction 3.5 a, (3 and k are arbitrary and 7 = — ^. In (3.2ii) we assume ci 7^ 
otherwise we get a classical reduction, and may set C2 = without loss of generality. Thus 
we obtain the following accelerating wave reduction 

u(x, t) = w(z) — eft 2 , z = x — \c\t 2 , 

where w(z) satisfies 

d 4 w dw d 3 w ( d 2 w 
W——7- + a— — —r + p 



Reduction 3.6 a and (3 are arbitrary, 7 = — ^ and k = 0. From (3.2iii) the following 
holds for arbitrary C2, and we may set c\ = without loss of generality. Thus we obtain 
the reduction 

u(x, t) = w(z)t — Cyt 4 , Z = X — |c2t 3 , 

where w(z) satisfies 

d A w dw d 3 w „ / d 2 w\ 2 d 3 w d 2 w f dw\ 2 dw 9 
w -TT + a -, — TT + P -nr ~ - w—^- - — + 4c 2 — + 12cS = 0. 



dz 4 dz dz 3 \ dz 2 

d 3 w d 2 w d 2 w ( dw\ 2 dw 



dz 4 dz dz 3 \ dz 2 J dz 3 dz 2 \ dz J dz 
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Reduction 3.7 k is arbitrary and a = (3 = 7 = 0. From (3.2iv) we get the following 
reduction 

U(x,t) = W(z) ± ± ^7=?/ 3 ± 7^2/ 2 + CQt, y = X±y/Kt, Z = X^y/Kt, 

where u;(z) satisfies 
d 4 u; 

^ dF 3 = 
This gives us the exact solution 

. . C3 4 o 2 C3 4 C2 q Cl 2 

u{x,t) = ^—=z +c±z- 3 + c 5 z + c 6 z + c 7 ± —^y ±—=y i ±—=y + c t. 

Reduction 3.8 k is arbitrary, a = — |, /3 = 2 and 7 = 0. In (3.2v) we may set ci = 
without loss of generality. Thus we obtain the following reduction 

u(x,t) = w(z)x 4 — \kx 2 , z = log(x) =F 

where w(z) satisfies 

d A w du; d 3 w / d 2 u>\ 2 d 3 w dw d 2 w 

4w—r - 6- r^r + 8 -r^r + 16^—^- + 58- 



dz 4 dz dz 3 \ dz 2 J dz 3 dz dz 2 

d 2 w d 2 w ( dw\ 2 dw _ 9 

+116w— 5- - k— T + 236 — + 776u>— + 672w 2 = 0. 
dz^ dz^ \ dz / dz 

Reduction 3.9 a = /3 = 7 = K = 0. From (3.2vi) and from (3.2vii) (04 = without loss 
of generality) we get the following reductions 

u(x, t) = w(x) + (c^x 3 + C2X 2 + C\X + Co) t 

and 

u(x, t) = w(x)t — (c^x 3 + C2X 2 + C\X + Co) 
respectively. In both cases w(x) satisfies 
d 4 w 

1? = °' 

These reductions have a common exact solution, namely 

u(x,t) = P 3 (x)t + Q 3 (x), 
where P3 and Q 3 are any third order polynomials in x. 

3.2.2 £ 2 + u = 0, not both a = 4 and 2£(f> u + £ u< /> = 0. 

The determining equations quickly lead us to require that both a = 4 and 2£4> u + £ u (fi = 0, 
which is a contradiction. 
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3.2.3 £ 2 + u = 0, a = 4, 2£<j> u + £ u = and /3 ^ 3. 

The determining equations give us that 7 = — |, k = and = 0. The invariant surface 
condition is then 

±i\/uu x + ut = 

which may be solved implicitly to yield the solution 

u(x, t) = w(z), Z = X^ly/ut. 

However, substituting into our original equation gives — — = 0, i.e. u(x, t) is a constant. 

az 

3.2.4 £ 2 + u = 0, = H(x,t)u- 1 / i , a = 4 and = 3, not all of # = 0, as = 0, 
1 + 2 7 = 0. 

For the determining equations to be satisfied, each of H = 0, k = and 1 + 27 = 0, which 
is in contradiction to our assumption. 

3.2.5 £ 2 + u = 0, = 0, a = 4, /? = 3, 7 = -| and k = 0. 

Under these conditions equation (^7|) which we apply the classical method to is identically 
zero. Therefore any solution of the invariant surface condition is also a solution of (|T|). 
Hence we get the following reduction 

Reduction 3.10 a = 4, (3 = 3, 7 = — | and k = 0. The invariant surface condition is 

±i\/uu x + ut = 
which may be solved implicitly to yield 

u(x, t) = w(z), z = x^iy/ut, 
where w{z) is arbitrary. 

4 Nonclassical symmetries (r = 0) 

In the canonical case of the nonclassical method when r = we set £ = 1 without loss 
of generality. We proceed by eliminating u x , u xx , u xxx , u xxxx and u xx tt in @) using the 
invariant surface condition (jig ) which yields 

U tt - ^(0^ + 00„) - 2j(u(f) x + -U00U + 2 ) - u(0 xa:a; + (j) u (f) xx + 2 02; + 00^ 

+40 U 2 MU + 50„00 a;u + 300,^02 + 3 uuu + 30 2 xuu + 300^ + ?><j)xu<t>x) 
/^(00M«^tt "I - 4"ruuu 

11% + 24>(j)t uu ut + 4>4>ttu + 2 1 «« + 30 n nu u 2 + 4<j) u (j)tuUt (32) 

+0n0« + 0z«^tt + 0x««^t + 2<f>xtuU t + 2(pt<PuuUt + ^.<t>t<t>tu + ^xtt) 
-a<p((f) xx + u 0a; + 00 2 + 2 un + 200j. u ) - (3((f> x + 00 M ) 2 = 

which involves the infinitesimal that is to be determined. As in the r / case we 
apply the classical Lie algorithm to this equation using the second prolongation pr^ 2 V 
and eliminate u# using (]32]). Similar to the nonclassical method in the generic case 
when 7^ the coefficient of the highest derivative term, uu is not necessarily zero, thus 
singular cases are induced. As in the previous section we consider the cases (i) [i = and 
(ii) \i 7^ separately. 
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4.1 Case (i) fi = 0. 

Generating the determining equations, again using symmgrp .max, yields three equations, 
the first two being 4> uu = 0, <fitu = 0. Hence we look for solutions like <fi = A(x)u + B(x, t) 
in the third. Taking coefficients of powers of u to be zero yields a system of three equations 
in A and B. 

aAA xxx + 2(3A 2 A XX + A xxxx + aA x A xx + 5f3AA 2 x + (in, 1,1; + WA XX A X 

+2>yA xx + 5AA XXX + aA 5 + WA 2 A XX + lOA 3 ^ + ,4 s + /M 5 + 15AA* (33) 
+47 A 3 + 2f3A x A xx + lOyAAa, + 4a,4 2 ,4^ + 6(3A 3 A X + 7aA 3 A x = 0, 

5a5A 2 + 2/M 2 5 :ra; + aBA xxx + 2/35A 4 + 135,4,4^ + 2aA 3 B x + KfyBAz 
+7^4^^ + 2(3A X B XX + a^S^j; + 2kAA x + a^ 2 5 xx + 15BA 2 A X 
+aB x A xx + 6~fAB x + 8jBA 2 + 2(3 A 3 B x + 2$B X A XX + 2aBA 4 + aAfl ssx 
• /W, + 6A BX B X + 2- /i.,,, + 5BA XXX + AB XXX + k,1.,,,. + 2/>\l' + A 2 B XX 
+A 3 B X + AA X B XX + 11BA 2 + ±(3BA 2 X + 6/3,4£ x ,4 x + 70^5^ + 7aBAA a:a: 
+2/?5AA x;c + W(3BA 2 A X + \2aBA 2 A x = 0, 

f3AB 2 x + 4 7j B 2 ^ + 2pB x B xx + 5B 2 AA X + a,45 2 + 6jBB x + 2kBA x 

+BAB XX + aB 2 A 3 + /?5 2 yl 3 + aB x B xx + 3aB 2 A xx + aBB xxx + BA 2 B X 
+3BA X B X - Btt + BB XXX + kB xx + B 2 A 3 + 3B 2 A XX + 2/?5A5 xa; (35) 
+aBAB xx + APBA X B X + 5aBA x B x + 2f3BA 2 B x + 2aBA 2 B x 
+4(3B 2 AA X + 5aB 2 AA x = 0, 

We try to solve this system using the dif f grob2 package interactively, however the expres- 
sion swell is too great to obtain meaningful output. Thus we proceed by making ansatze 
on the form of A(x), solve ( |3~5"1) (a linear equation in B(x,t)) then finally J34l ) gives the 
full picture. Many solutions have been found as ( |33[ ) lends itself to many ansatze through 
choices of parameter values. We present some in § 4.3. 



4.2 Case (ii) // ^ 0. 

The nonclassical method, when the coefficient of uu is non-zero, generates a system of 
three determining equations. However, far from being single-term equations the first two 
contain 41 and 57 terms respectively, and the third 329. The intractability of finding 
all solutions is obvious. To find some, we return to our previous case and look for <f> = 

,3^j35|) which we tackle in 
the same vein as previously. Some solutions are presented in § 4.3. 

As mentioned in the start of this section, singular solutions may exist, when the 
coefficient of uu equals zero, i.e. when 

1 - <f><t>uu - 4>u ~ <t>ux = 0. 

This may be integrated with respect to u to give 



A(x)u + B(x,t). Three equations then remain, similar to ( 



u - 4>4> u - 4> x = H(x,t). 
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If 4> satisfies ( [4.2|) then the coefficients of u 2 and u t in (|32]) are both zero. Since no 
it-derivatives now exist in (^2|) what is left must also be zero, i.e. 

(2 7 + a)<f - aH x (j> + (2 7 + + l)u 2 

+ [k- (2 7 + 2/3 + 1)# - u + I3H 2 -kH- H tt = 0. 

Thus we need to solve ( |4.2[ ) and Q4.2j ). Note that once we have found cj)(x,t,u), the 
related exact solution is given by solving the invariant surface condition, with no further 
restrictions on the solution. The following are distinct from each other and from solutions 
in § 4.3. 

Case (a) 7 = — ^, a = 1 and (3 = k = 0. In this case 4>(x,t,u) is given by the relation 

u - <j>4> u - 4> x = c\t + c 2 . 

For instance, if (j>(x, t, u) is linear in u we have the exact solution 
u(x, t) = w(t) coshfx + A(t)} + B(t) sinh[x + A(t)] + c\t + c 2 , 
where w(t), A(t) and B(t) are arbitrary functions. 

Case (b) a = —2<y, (5 = — 1 — 27 and 7 / — |. In this case </>(x,t,u) is given by the 
relation 

u - 4>4> u - <p x = -— —. 

I + 27 

For instance, if (f>(x, t, u) is linear in u we have the exact solution 



u(x, t) = w(t) cosh[x + A(t)j + B(t) sinhfx + A(t)] 



K 



l + 2 7 ' 

where w(t), A(t) and B(t) are arbitrary functions. 
Case (c) a = —27 and j3 = — 1 — 27. In this case 

^(x, i, u) = ±u — KX =F (— ^Kt 2 + Clt + C2) ± K, 

and so 

t) = w(t) exp(±x) ± kx + (—^Kt 2 + c\t + c 2 ) , 

where w(t) is an arbitrary function. 

Case (d) 7 = — |, a = 1 and = 0. In this case 

Kit - i^li - kH - H tt 
<p{x,t,u) = , 

where H x (x,t) 7^ and also H(x,t) satisfies the system 

kH xx + H 2 -k 2 = (k 2 - Hi) H u ~ 2 K H 2 t + k 2 (k 2 - Hi) = 0. 

We have assumed that k 2 — Hi 7^ 0, for a different solution to (c). This yields 
u(x, t) = [w(t) — 2kx] sinh z cosh z 

+ cosh 2 z [4k log(cosh z) - n 2 t 2 + 2c 3 t + 2c 4 - 2k + 2cf\ - cj, 

where z = \{x + c\t + c 2 ) and w(t) is an arbitrary function. 
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Case (e) (3 = — 1 — a — 47, 7 7^ I and a + 27 / 0. In this case 



(j)(x, t, u) = ± ( u + 



1 + 27 / ' 



and so 



t) = w(t) exp(±x) — 



l + 2 7 ' 

where w(t) is an arbitrary function. 

Case (f) 7 = 0, a = —2 and (3=1. In this case (p(x,t,u) satisfies 
-2(f) 2 + 2cpH x + 2u 2 + 2ku - 2Hu + k 2 + H 2 = 

and H(x,t) satisfies the system 

H xx + k + H = 0, H tt + k 2 + kH = 0. 

Then 

u(x,f) = i(A 2 + J B 2 ) 1/2 sinh[±x + «;(<)] - « + ^(Asinx + Bcosx), 
where A(t) and B(t) satisfy 



d 2 A 



+ kA = 0, 



d 2 B 



+ K.B = 



dt 2 ' dt 2 

and w(t) is an arbitrary function. 

Case (g) 7 = 0, (3 = — 1 — a, k = and a = (ci — l) 2 /ci, where c\ / 0, 1. In this case 

it 1 \ . c 2 t + c 3 . . 
<p(x, r, uj = u H — — exp(— c\x), 



and so 



u(x, t) 



ci-1 



ioi(i)e x - i(c 2 t + c 3 )xe 2 



if a = -1, 



t«2(t)e ie + (1 - c 2 ) i (c 2 t + c 3 )e- cia; , if ci^-l, 



where and u> 2 (£) are arbitrary functions. 

Case (Ti,) 7 = 0, /? = -1 - a, cf + a 2 ci + 2c\ + 4aci + 1 = 0, a / 0, -2 and ci / 0, 1. In 
this case 0(x, t, u) satisfies 

ac/) 2 - aH x (j) -au 2 + u[k + H(l + 2a) - H xx ] - kH - (a + 1)# 2 - fl« = 0, 

where H(x, t) satisfies the system 

H tt + Ci kH = (a + 2)fl a . ± (1 - ci)(ff + k) = 0. 

Then 



u(x, t) 



1 + 2a + a 
2a 



w(t) +g(x,t) 



exp 



(ci - l)x 
a + 2 



- K, 
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where w(t) satisfies 
d 2 w 

and g(x,t) satisfies 

2 2a(c\ — 1) 2 

a 9 x H "Ti — - «(ci + 1)5 + (aci + ci + l)w{t)g x 

a + 2 

- Cl (l + « + ^ + Cl(Cl " 1)(a + 2) ^ 2 (t) = 0. 

4a 

4.3 Exact solutions 

In this subsection some exact solutions are presented. The infinitesimal <f>(x,t,u) is given, 
possibly up to satisfying some equations, and then the solution, found by solving the 
invariant surface condition (jl6|). 

4.3.1 7 = and <j> = ~ + Hi(t)x + 3H 2 {t)x 3 + + H A {t)x 2 - a . 

Solving the invariant surface condition gives 

u(x, t) = < 



xw(t) + H x {t)x 2 + H 2 (t)x 4 - H 3 (t) + g4 (t)xS " , if a/ 2, 



2 - a 

2 , u /.a„,4 



^(i) + + H 2 (t)x * - H 3 (t) + # 4 (f)x log x, if a = 2. 

Various types of solution are found, as seen in Table 4.1. The Hi{t) are obtained from by 
the determining equations, w(t) by substituting back into (||). 

4.3.2 <p = B(x,t). 

Case (a) 7 = 0. In this case B(x,t) = 4Hi(t)x 3 + 3H 2 (t)x 2 + 2H 3 (t)x + ffii(i), where 
H2(t), ff 3 (i) and H 4 (t) satisfy 



d 2 # 



di 2 

2 



24(6/3 + 4a + l)Hf = 0, 



d 2 K 



24(6/3 + 4a + l)HtH 2 = 0, 
<l 11:1 24(2/3 + 2a + l)#i# 3 = 18(2/3 + a)^ 2 + 12fsfli + 12/x%^, 



di 2 
f 3 



dt 2 " M ' " ' ; 1 d ^ H ' ^ 2 ' 1 ' ^ dt 2 

j2rr j2ij 

' ' 1 24(a + l)H x H A = 12(2/3 + a)H 2 H 3 + 6kH 2 + 6fi- " ' 



Then 

u(z,t) =w(t) + Hx{t)x A + H 2 (t)x 3 + H z (t)x 2 + H A {t)x, 
where uu(t) satisfies 

^ - 2AH lW = 2kH 3 + 6aH 2 H 4 + 4/3tff + 2^^^. 
di^ d^ 
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Table 4.1 



Parameters 


i?i(t) and satisfy 


u 2 , /J g 


^ = 4k, H 2 = H 3 = ^ = 
d 2 w „ TT 9 

32 ^ = -m\ 


a = ^, k = 


H x =f iH 2 = H 3 = 
^-72(1 + 2/3)^1 = 

J 2 IT 

16 , 4 3(480/3 + 303)H 2 H 4 = 
di^ 

8^ - 288H 2 w = (50/3 + 5)tff 


9 

3 — a 
k = 0, a / 3 


ffl = ^#2 = #3 = 

(a - 3)^9^ + 24(2a + 3) (a + 1)F 2 2 = 

d 2 4 +3(a + 2)(a + l)(a 2 a A)H 2 H A = 
d 2 w 

— T - 24(q + l)H 2 w = 
dt z 

^ - 72H 2 w = 90H 2 H 4 (if a = 2) 
dt z 


a = -2, /3=§ 


^ = §«, ^ 2 = ^ = 

d 2 rYo 

5 d(2 -24ff|=0 

5—^ - 120H 2 H 3 = -25k 2 
dt z 

^ + 2AH 2 w = -30H 3 H 4 



Case (b). In this case B(x,t) = Hi(t) + 2H 2 (t)x where H\{t) and H 2 (t) satisfy 



d 2 H, 
dt 2 

d 2 H, 
dt 2 



9 o 

12 7 H 2 



0, 



12jH 2 H 1 = 0. 



(36) 
(37) 
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Then 

u(x,t) =w{t) + Hi(t)x + H 2 (t)x 2 , 
where w(t) satisfies 
d 2 w 



dt 2 



- 4jH 2 w = 2kH 2 + 4(6 7 /i + 0)H$ + 2-yHf. (38) 



Case (c) (3 = 1-a. In this case B(x, t) = cHi(t)e cx + cH 2 (t)e- cx + H 4 {t), with c 2 = -2 7 
and where Hi(t), H 2 (t), H 3 (t) and H 4 (t) satisfy 

^ = 0. (39) 
d 2 H 

(1 + 2 7 /x) - 2 7 c(2 - a)H A H! + 2 7 k#i = 0, (40) 
d 2 H 

(1 + 2 7A i)— + 2 7 c(2 - a)F 4 iJ 2 + 2 7 kF 2 = 0. (41) 

Then 

u(x,t) = w(t) + H^ty* - H 2 (t)e~ cx + H 4 (t)x, 
where w(t) satisfies 
d 2 

-^ = 2 7 H 2 -W 7 2 (l-a)H 1 H 2 . (42) 

Case (d) a = 2 and (3 = -1. In this case B(rr,t) = cH^t)^ + cH 2 {t)e~ cx + 2H 3 {t)x + 
i? 4 (t), with c 2 = -2 7 and where ffi(t), # 2 (i), #3(0 and # 4 (i) satisfy 



d 2 H : 
dt 2 

d 2 H 4 
dt 2 



'' 12 7 #f = 0, (43) 
12 7 H 3 H 4 = 0, 



d 2 i? 

(1 + 2^)-^ " 12 7#3#i + 2K 7 i7x = 0, 
d 2 H 

(l + 2 7 /i)-^-12 7 # 3 # 2 + 2K 7 tf 2 = 0. (44) 

Then 

u(x,t) = w{t) + H^Qe * - H 2 (t)e~ cx + # 3 (t)x 2 + H 4 (t)x, 
where w(t) satisfies 

^ - 4 7 H 3 w = 2kH 3 + 2 7 #J + 4(6 7/ u - 1)F 3 2 + 16 7 2 #i# 2 . (45) 



90 



P.A. Clarkson and T.J. Priestley 



4.3.3 7 = f3 + a + 1 = and <f> = R(u + H^t) + H 2 (t)e Rx + H 3 (t)e m + X + tf 4 (t)e m - x ). 

Here m± = — ^i?(2 + a db n), n = a (a + 4), with R ^ ±1 a non-zero constant. Solving 
the invariant surface condition yields 

u(x,t) = w(t)e Rx - Hi(t) + RH 2 (t)e Rx 

_2H 3 (t)_ (_i Rx , 2 + + xi _ 2ff 4 (t) (i Rx <2 + a _ n \\ 
4 + a + n ;J 4 + a-n L 2 J 

The solutions are represented in Table 4.2 

The equations that the various H{(t) satisfy in this subsection are all solvable, and 
the order in which a list of equations should be solved is from the top down. The only 
nonlinear equations all have either polynomial solutions (sometimes only in special cases 
of the parameters) or are equivalent to the Weierstrass elliptic function equation (^). The 
homogeneous part of any linear equation is either of Euler-type, is equivalent to the Airy 
equation [||, 

or is equivalent to the Lame equation p3] 

^-(t)-{k + n(n + l)p(t)}H(t) = 0. (46) 

The particular integral of any non-homogeneous linear equation may always be found, up 
to quadratures, using the method of variation of parameters. 

For instance consider the solution of 4.3.2 case (b) above. There are essentially two 
separate cases to consider, either (i) 7 = or (ii) 7 / 0. 

Case (i) 7 = 0. The functions H\(t) and H2(t) are trivially found from (^) and ( |37l) to 
be Hi(t) = c\t + c 2 and H 2 (t) = c 3 t + c 4 , then (|38|) becomes 

d 2 w 



2 K (c 3 t + c 4 ) + 4/?(c 3 t + c 4 
ated twice to yield 

Q 

Jc 3 t + c 4 ) 3 + —2 (c 3 t + c 4 



dt 2 

which may be integrated twice to yield the exact solution 



u(x, t) 



ft 1 \ 4 

„— (c 3 £ + c 4 ) 

>rr i 3c 3 if c 3 /0, 

+c 5 t + c 6 + (c\t + c 2 )x + (c 3 t + c 4 )x 2 , 

. (kc 4 + 2/3cl) t 2 + c 5 t + c 6 + (cit + c 2 )x + C4X 2 , if c 3 = 0. 



Case (11)7^0. Equation ( |36| ) may be transformed into the Weierstrass elliptic function 
equation (p9|), hence H 2 (t) has solution H 2 (t) = p(t + to; 0, g 3 )/(2j). Now Hi(t) satisfies 
the Lame equation 

d 2 H 

-^- L -6p(t + t ; 0,g 3 )Hi = 0, 
which has general solution 

rt+t ds 



Hi(t) = cxp{t + t ;0,g 3 ) + c 2 p(t + 1 ; 0, g 3 ) J 



P 2 {s;0,g 3 ) 
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where c\ and c 2 are arbitrary constants. Now w(t) satisfies the inhomogeneous Lame 
equation 

-^-2p(t + t ;0,g 3 )w = Q(t), (47) 

where Q{t) = 2KH 2 (t) + 4(67/x + (3)H 2 {t) + 2jHf(t), with H x (t) and H 2 (t) as above. The 
general solution of the homogeneous part of this Lame equation is given by 

W CF (t) = C 3 Wi(t + t ) + c 4 w 2 (t + t ), 

where C3 and C4 are arbitrary constants, 

wi(t) = expj-tC(a)} 7-—, w 2 (t) = exp{tC(a)j- 



a(t) ' ^ — — J a(t) 

in which C(z) and cr(z) are the Weierstrass zeta and sigma functions defined by the 
differential equations 

^ = -pO), ^ log cj(z) = CO) 
together with the conditions 

lim ( C(z) --^1=0, lim ( ^ 

z^O \ Z J z-+0 \ Z 

respectively (cf. [f7(]]), and a is any solution of the transcendental equation 

Pifl) = 



i.e., a is a zero of the Weierstrass elliptic function (cf. |43| |, p. 379). Hence the general 
solution of (E^) is given by 



w(t) = C 3 Wi(t + t ) + C 4 W 2 (t + to) 



1 



ft+to (48) 

+ \ / [wi(s)w 2 (t + to) - wi(t + t )w 2 (s)] Q(s) ds, 
W{a) J 

where W(a) is the non-zero Wronskian 

W{a) = w\w' 2 - w[w 2 = -a 2 (a)p'(a) (49) 



and Q(t) is defined above. We remark that in order to verify that ( [48| , [49| ) is a solution 
of (|47| ) one uses the following addition theorems for Weierstrass elliptic, zeta and sigma 
functions 



C(*±t) = C(*)±C(t) + 5 



P'(s) =F P'(t) 



P(») - P(*) 
(1(5 + t)a(s - t) = -a 2 {s)a 2 {t)[p{s) - p(t)] 

(cf. m, p.451). 
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Table 4.2 



Parameters 


Hi(t) and w(t) satisfy 


a = -4 
/3 = 3 


h h d2Hl - n 

H a -H A - ^ -0 

(l - fiR 2 )-TT^ + R 2 (R 2 Hi - K)H 2 = 
at z 

(1 fiR 2 ) d *™ + R 2 (R 2 H 1 k)w = 2kR 2 H 2 AR^H^ + liiR 2 ^** 2 
at z dt z 


a arbitrary 

J 2 31 2 
6 2 12 


h h d2Hl n 

2 ~ j ~ H/2 ~° 

(4 - n(2 + q ± n) 2 R 2 )-^- + i? 2 tf; 
x [i? 2 F!((Q 2 + 4a + 2)(2 + a ± nf - 4) - «(2 + a ± n) 2 ] = 

(l vR 2 ) d d ™ + R 2 (R 2 H 1 k)w = 


a = -3 

/3 = 2 


H2 ~ dt 2 "° 

J 2 TT 

(2 + ^ 2 (1 + iV3))-^ - HiHsR 4 ^ - iV3) + K# 3 i? 2 (l + iv/3) = 

J2 T_r 

(2 + /xfl 2 (l i>/3)) 2 4 H 1 H A R A {l + \V3) + KHiR 2 {l i>/3) = 
(1 M^ 2 )^™ + R 2 (R 2 H 1 k)w = 6H 3 H A R 4 


a = — 1 

(3 = 


H2 ~ dt* "° 

(2 + ^ 2 (1 - 1^3))^^ - H i R z R \ l + i>/3) + /tF 3 i? 2 (l - iV3) = 

J 2 rr 

(2 + /xi? 2 (l + h/3)) 2 4 ^iF 4 i? 4 (l iV3) + KF 4 i? 2 (l + i\/3) = 
(1 nR 2 ) d d ™ + R 2 (R 2 H! k)w = 
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5 Discussion 

This paper has seen a classification of symmetry reductions of the nonlinear fourth or- 
der partial differential equation (|]) using the classical Lie method and the nonclassical 
method due to Bluman and Cole. The presence of arbitrary parameters in (||) has led 
to a large variety of reductions using both symmetry methods for various combinations 
of these parameters. The use of the MAPLE package diffgrob2 was crucial in this 
classification procedure. In the classical case it identified the special values of the param- 
eters for which additional symmetries might occur. In the generic nonclassical case the 
flexibility of diffgrob2 allowed the fully nonlinear determining equations to be solved 
completely, whilst in the so-called r = case it allowed the salvage of many reductions 
from a somewhat intractable calculation. 

An interesting aspect of the results in this paper is that the class of reductions given 
by the nonclassical method, which are not obtainable using the classical Lie method, were 
much more plentiful and richer than the analogous results for the generalized Camassa- 
Holm equation (|2|) given in [24\. 

An interesting problem this paper throws open is whether (|l|) is integrable, or perhaps 
more realistically for which values of the parameters is (|]) integrable. Effectively, in 
finding the symmetry reductions of (|l]), we have provided a first step in using the Painleve 
ODE test for integrability due to Ablowitz, Ramani and Segur g]. However the presence 
of so many reductions makes this a lengthy task and so the PDE test due to Weiss, Tabor 
and Carnevale |6^] is a more inviting prospect. It is likely though that extensions of this 
test, namely "weak Painleve analysis" |5ll] and "perturbative Painleve analysis" p£ 



will be necessary (for instance see [flOf ). We shall not pursue this further here. 

The FFCH equation (|ll]) may be thought of as an integrable generalization of the 
Korteweg-de Vries equation @. Analogous integrable generalizations of the modified 
Korteweg-de Vries equation 

qt = q x + qxxx + 3iq 2 qx, 

the nonlinear Schrodinger equation 

m + qxx + \q\ 2 q = o 

and Sine-Gordon equation 

q xt = sin q 
are given by 

ut + vu xxt = u x + q xxx + 7[(n + vu xx ){u 2 + vu x )] x , (50) 
iut + iu x + fiUxt + u xx + ku\u\ 2 — iKfi\u\ 2 u x = 0, (51) 
u xt = sin(n + fiu xx ) = 0, (52) 
respectively, where /i and k are arbitrary constants |32, |33|. 



Recently Clarkson, Gordoa and Pickering 17] derived 2 + 1-dimensional generalization 
of the FFCH equation (O) given by 



r ^"iiU xxxx ~\~ ^xy^xxx & {^2^ 3 y r ^ , xx ~t~ ^x^xy) ^xxxt Q-^xt — 0j (^^) 



2 ""[I " r.r.r.r 
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where a is an arbitrary constant. The FFCH equation (11) and is obtained from (53) 
under the reduction d y = d x , with v = u x . The 2 + 1-dimensional FFCH equation (11) 
has the non-isospectral Lax pair 

4ip xx = [a - X (u xxx - au x )] ip, 

Ipt = A" 1 1py ~ \Uyi) X + \u X y1p, 

with A satisfying A^ = AA^. Clarkson, Gordoa and Pickering |17j also derived a 2- 
component generalisation of the FFCH equation (|ll|) in 2 + 1-dimensions given by 



t^xxxt OLU x t — 2^"y'^' xxxx / U / xyU xxx ~\~ OL [^^Uy'U'xx ~\~ U x U X y) KU XXX y + Vy, 

V t = ~VU X y ~ \v x Uy, 

which has the Lax pair 

4(1 + K,X)ip xx = [a - X (u xxx - au x ) - X 2 v\ ip, 
ipt = A"Vy - \u y %l) x + \u xy tp, 

where the spectral parameter A satisfies X y = XXt- 



(54) 



We believe that a study of symmetry reductions of (5C ,51 ,52, 52 : 54) would be interesting, 
though we shall not pursue this further here. 
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Appendix A 

In this appendix we list the determining equations that are generated in § 3.2 in the generic 
case when £ 2 + u ^ 0. 

tu = 

a£ 2 <pu - 2a£& + A<\> uu u 2 - acj) + 4£ 4 uu + a<f) u u + 8{, 2 (p uu u - 2a£ 2 £ x = 
f3t 2 <p u + 3(j) uu u 2 + 3^ uu + Ge^uuu -(3(1) + P<t> u ii - 2/3£% - = 
I2£, 2 <fr uuu u + 3acf> uu u + 6<fi uuu u 2 + 2(3(f> uu u + 2f3£ 2 4> uu + 3a(, 2 (j) uu + 6£ 4 </w = 
7£ 2 £xxU - 10£gu - 2i 2 i t (\) u - A^ xt u - ai 2 (p x - 2£<j><j> u - - ac\> x u + 2&0 u u 

+2£ 3 # un + 5& x <f) + 2£(j) tu u - &i 2 <j) xu u + 4£%£ x - 2£& 2 + 6£ xx u 2 + 2( 3 <p tu + t tt u 
H A ixx + i 2 ttt + 4££,x<t>uu - 4(p X uU 2 + 2i^ uu u - 2^4> xu - 4C 3 ^t - it<P = 

(fruuu +~ CXtfiuuuU ~i~ 2£ (fruuuuU ~i~ C (fruuuu ~i~ 4'uuuu'U' — 

3a£ xx u - A(3£ 2 (() xu - 12(p xuu u 2 - 3acf) xu u + 6£cf) tuu u + 3a£, 2 £ X x - 18£ 2 4> xuu u 
+^it(puuu - Af3(j) xu u - 6£(j)<j)uu + 2(3£ xx u + Q^4>4>uuu - ^^ix<t>uu - 3£ 2 £t4>uu 
-3a^ 2 (j) xu + 12£ 3 4> u (j) uu - 3^ x 4> uu u + 6^(f)(t>uuuU + 12£<f> u <j) uu u + 2(3£ 2 Z, XX 
-6£, 4 4>xuu + 6£, 3 4>tuu = 
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2£, 3 4> X u4>xx ~ Z, Z £xx4>xx + 2£&# xxu - 2££ t (j) t + £ xx 4> xt U + ^4> xxu 4>x + 2£xxt<t>xU 

+4> 2 (f) U uU - 2^ t 4> x u + 2£ s 2 4>4> tu - £ xx # x - 4>t4>xxuU - 2£ x <j) xxt u - 2££ t £ xx (f> x 
+2(fxf) tu u - K(f> xx u - 2<j) x t<j) 2<f>4>xxtuU + 4£€t<f>xu<l>x + 2££ t <t>xxt 

-2£ 2 4> xt (pxu + ^x4>tu - 2£ 2 4>4> xxtu - 4£ 2 4> xtu (t) x + ^i 2 i xt 4> xx - i 2 $ 2 4>xxuu 

+<t*l>u<t>xx - (?K<j>xx + £ 2 txx4>xt - £ 2 <Pt4>xxu + £(f>(t>x - 24>(/) 2 xu U - 2^ 2 (j) uu 4) 2 x 

-2£ x <jxj) xx + 2£ t (f) xxx u + 2t > 2 £ tX 4> xu (j) x + £ 2 4> 2 (f) uu - 4> 2 4> xxuu u + 2<jxf> xu <f) x 
-i 2 4>xxx X u + 2^ 2 $ >xxt (f) x - 2£ 3 £ r </> x - 2^ 2 4>4> 2 XU - 2(p tu (j) xx u + ^ 2 cj) xx u - 2^<f> xx u 2 

+ 2£ 2 £ x <f>t - 2t, 2 (j)tu<t>xx - 4>4>u(t>xxuU - £4>4>xxx - 2^ t <f>4>u - 2<f> uu <f>lu - 4<f> xtu (f) x U 
-2(j>(f>uu4>xxU - 4(f)(t)xuu4>xU - i 2 4>xxtt + i 2 4>tt - 4>xxxxU 2 + 2££ x (f) xxx U + £, xx (f)4>xuU 
+£, 2 £xx<t>u4>x ^x^xxfixU ^CxxtpxxU ~i~ ^ < fixxu4 > x'U + £xx4>u<PxU ~ 2£ x (f)<t> xxu U 

-2^4>^uu4>xx + 4£ x (j)(f> u u + 2£4> xu <f> xx u - 2£ x 4> u (j> xx u - 2{ t 2 ^(j) xx u - 2£ x (p xu (f> x u 

-Z, 2 £,x£,xx4>x + 2i 2 i x ^ u + 2££ t (j) u (/) xx - 4££t£x<f>xx - ^ 2 <t>4>xuu(f>x - (pxxttU 

+4> tt u - 4> 2 4> u + 4>4>xxt - 4><f>t + 2 4>xxu - 2i 2 ^ xu ^ x - i 2 ^<\) u ^ xxu - 4££ x x u 

+i 2 ixx<P4>xu - 2(p u (j) xu <p x U = 
2£ 3 £xx(pu - 4£ 2 £xxxU - 8£&£ 2 - 3Z 3 Zxtxx + 4£€ttxt + 4> 2 <t>uuuU ~K(j) + 4£ xt £ x u 
+2cf) tu (f) u u - 2£&u - ^cj) xtu u - ^Hxuu + 2[3£ 2 <j) xx - 4£ 2 £ x « + 4£ x7 u 2 + 2£ 2 7 
-^xt(t> u u + 2(j)4> tuu u - 4£ 2 £ x (t>tu - ±i 2 ixt4>u - ^<i>uu4>x - S&^uU + 2£ x 4> 2 u u 
-2£&k - A^4> u <j) xu + 7£t€xxU + i 2 4>t4>uu - 4£€xx</> + <Pt<Pu U u + t 2 4> 2 <Puuu 
+2££ xxt u - i 2 itixx + 2i 2 it4>xu - 8£t<f>xuU + 6£ 3 £ x xtt - 2^ x 4> tu u - 2££ t <f><f> uu 

+2^ 2 4>tu4> u - 2^t4>1 + 7£ 2 <t>xxu,u + %££,t£x<t>u - ^4>u4>x U u - 2£ x <jxj) uu u - \ii x 4> xu u 

+2j34> xx u + 2i 2 4>4> tuu + 2H xx 4> u u + fc/) ttu + 8gu - 4f(/) xtu + 2£ 3 £ xxi + 5^ xu 
+2£ xt <f> - 4££ t ju + i A 4> xxu + 5££ x tixxU - 5£,4>uu<f>xU - ^4>4>xuuU - 2££ t <j) tu + 4> Uu u 

-4> 2 <t>uu + 3^ 2 00 U U „ + 3<f)(f) u (f) uu U - (fxfrl - <f)<f>tu + 2£ X KU - ^i x <\>4>uu + ^x4>4>u 

+64>xx U u 2 - 4£ xxx u 2 - 2^ 2 ^ x « - 2i A i x - 2i xtt u - Ai 2 x d) + £ 2 </> + 8^ 2 ^ a = 
2i4>tuuuU - h£ z £ x 4>uuu - 3a(j) xuu u + 2£, 3 <jxj) uuuu + 6^ 2 M n + 3(, t 4>uuuU + 6£ 3 4>u4>uuu 
+6£, 3 <Plu + 2£?(f)tuuu - 2^4>4> uuu - f?it§uuu - 3a£ 2 4> xuu + 6£,4> u (j) uuu u 

-\-2^(j)(f) uuuu 'U 6£ (frxuuuU £,£,x<t'uuu'U' ^(ftxuuuU 2£ 4>xuuu = 

^ 2 <Ptu4>uu + £, 2 4>t4>uuu - a£, 2 £ xxx + 4:£ 2 <j) tuu <j) u - 6£ xt (f>uuU - 270 + 44> tu 4>uuU 

+ 2(fxf> tuuu ii - a£ xxx u + 5^(f)(f) xuu + 2£ 2 £ t (f) xuu + 4£ x cf)(/) uu 
-<±i 2 ix$tuu + 2£ 2 7^ n + 4^ 2 ^ 2 ^ UM - 5^ 3 (j) uuu 4> x + 27^ MU n 2 + £ 2 K<p uu + 6£ 3 £ x <p xuu 

-8£ 3 (pu4>xuu ~ ^ 3 4>(t>xuuu - 2^ t <f>tuu + 4<f>tuu4>uU + <f>t<f>uuuU ~ 14:£ 3 <f> xu (f) uu 
+ 2^ 2 (f>(t>tuuu - Ht<t>xuv,U - §i 2 £,xt4>uu + ^ 3 ixx4>uu + 3a<f> xxu u - 2£ x (j) tuu U 
-£ 2 4> uu u + ^4>4>uu u + k4> uu u + 4£ r7 « + 4:(pl<j)uuU - 8£,4> u (p xuu u - 4:£<f) xtuu u 
+i 2 4 >2 4 ) uuuu + 7i 2 <t>xxu U u - 14^cj) xu (j) uu u - 2£ x c/) u 4> uu u + 2£ 2 j<f> uu u + 6££ xx (f> uu u 

+^ 2 <Pl4>uu - 4£&7 - <j) 2 <puuu - 4>4>tuu + 6<PxxuuU 2 - ^ 3 4>xtuu - 2£&# uuu 

-4:£<jxj) xuuu u - ^l4>uuU - 5£<j) uuu (j) x u - ^4> uu - 8£, 2 £ >x <j) u (f) uu + 5<f)<t) u <p uuu u 
u'Puuu + 8iitix4>uu — 2^ x 4>4> + 3at 2 4> xxu + 

xxuu 3<fi(f)u4>uu 

-4£€x<f>xuuU - 4£ 2 (, x 4>4> uuu + i 2 <\>ttuu + (pttuuU + 4^ 2 00 2 n - 6&t4>u<Puu = 
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tit (fixuU 

+ 2<j) t (f) xuu u + £, 2 £,t4>xxu + ^-l4>xU 

+MxxxuU 2 + 24>xttuU + §(, 2 <j)<i>u<t>xuu ~ ^ 2 ^x<t>4>xuu + 4£, 2 <p xtu 4>u + ^ 2( P4>xtuu 

-2(f) 2 (t) xuu + 2(, xx 4>(j) u + 2£ 2 2 ;EU - 4f(f> uu 4> xx + 20 2 xn u + - 2£ 2 £ t u 

+2£0 iu u + + 2£?(j)(j) uu - 2£ 3 <p u (j) xxu - 4<pcf) u <p xu - W£, 2 ^xt4>xu - ^ t <t>xtu 

+^ 2 0tuu(px + 2£ 2 6£r + ^(j) xu U 2 + H 2 ixttxx + a<PxxxU + ^ x (f> + 2£ 2 (j) 2 (p xuuu 
-4{,£t(t><t>xuu + €>i 2 <t>tu<t>xu + 5££xx4>xuU + 2K<fi xu U - W£, xt <f>xuU + ^ 2 -f4> x - £ XX KU 

-3ixixx4> + 2i 2 i xx £ x + 2£ 2 n<j) 

xu ~i~ £ Cx^xxu xt^uu ~i~ 2£ ^(frxu £,£,xxx4 > 

-2£ 2 £ xxt <f> u - 2(fxp uu <f> x - i 2 i xx (f) 2 u + 2£ t (j) u u + 2£ t £ xxx u + 5£, 3 £ xx 4>xu + 2££ t £xxt 

~i 2 i 2 xixx + ^4>4>uuu4>xU + §4>u4>uu4>xU - 2<f)<f> xtu + 8£, 2 (fxt) xu (f> uu - 2£ x £ xx (f> u U 

-?>i 2 ixx(t>tu - 2i4> xxtu u - a 2 xx u - i xxtt u - ?>i 2 i xx 4>4> uu + 2^^ - ^£ 2 XX 

-2(?(j)<j) xxuu - £ 2 £ XXXX U + 2£ 2 <f> t <f>xuu - ixx4>l,U + §ix4>4>xu - 2£ >xxt (j) u U 

-4£, 2 £ x 4> x tu - 2£00 u - £ 2 £ xx k + b£, xt t,xxu + a^4> xxx - 8^ t 4>u<f>xu - 4££ 2 w 

-^4> 2 xu + §4>4>u4>xuuU + 8<j)<p xu <j) uu U - ^i 2 i x ^u^>xu + ixxt4> - 2£?(f> xx tu 

-l^l<j) xu u + A(f) xtu <j) u u - 2££ t 2 - 4££ t cf) uu (f) x - 3£ xx (f)tuU + 24> 2 4> xuuu u 
+12££ t £x4>xu - ixxxxU 2 - 2£ x 4> uu 4> x u - 2£ 3 £ 2 + 2£ 2 <p xttu + 4£ 3 £ r u - 6((f> xuu 4> x u 
+2^4> uu u - 2£4> u <j) xxu u - 2^4>4> xxuu u + ^ 2 i4> xu u + 2^ 2 ^ xx <j) u - 7£t4>xxuU 

~i 2 ixxtt + ^ 2 4>xxxuU - 7££x4>xxuU + ^ 2 <t>(puuu4>x ~ ^Utixixx + £, 2 £,tt 

-^4>uu(pxxU - <oi 2 ix$uu$x + 6£ 2 4>u4>uu4>x ~ 2i 2 i xx ~ju - 3£, xx (j)<j) uu U + 2££ X £ XXX U 

-60 + ^x4>u<PxuU + 4£,£t£,xx4>u = 
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